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ABSTRACT 

We present a new and improved model independent amplitude analysis of reac- 
tions n~^n^ — > TT~^TT~p at 5.98 and 11.85 GeV/c and n~p^ — > TT~TT~^n at 17.2 GeV/c 
measured with transversely polarized targets at the CERN Proton Synchrotron. 
For dipion masses below 1000 MeV the pion production process is described by two 
5'- wave and six P-wave production amplitudes. Previous analyses suffered from the 
presence of unphysical solutions for moduli of amplitudes or cosines of their relative 
phases, causing uncertainties regarding the signal for scalar state / = 0"'""'"(750). 
To remove the unphysical solutions we use a Monte Carlo approach to amplitude 
analysis. In each (m, t) bin we randomly varied the input spin density matrix 
elements 30 000 times within their experimental errors and performed amplitude 
analysis for each selection. Unphysical solutions were rejected and the physical so- 
lutions produced a continuous range of values for moduli, cosines of relative phases 
and for partial wave intensities. A clear signal for cr(750) resonance emerges in all 
four solutions for S'-wave intensity Is at 5.98 and 11.85 GeV/c and in both solutions 
for S'-wave amplitude |5'pE at 17.2 GeV/c. Its tt+tt" decay width is estimated to 
be in the range 200-300 MeV. We find significant suppression of p° production in 
the amplitudes |iVpS and, at 17.2 GeV/c, in |LpS. The mass dependence 

of amplitudes \L\^T, and \L\^T, shows unexpected structures within the p° mass 
region which correlate the mass spectra corresponding to opposite nucleon spins. 
These features of P-wave moduli reveal the essential role of nucleon spin in pion 
production process and contradict the factorization hypothesis. Our results empha- 
size the need for a systematic study of pion production on the level of amplitudes 
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in a new generation of dedicated experiments with spin at the recently proposed 
high-intensity hadron facihties. 
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1. INTRODUCTION 

The pion production process ttN — > Tr+Tr^A?" has always served to develop our 
ideas on the dynamics of hadron collisions and hadron production. In 1978, Lutz and 
Rybicki showed^ that measurements of pion production in meson-nucleon scattering 
on transversely polarized target yield in a single experiment enough observables that 
almost complete and model independent amplitude analysis can be performed. In 
the kinematic region with dimeson masses below 1000 MeV the dimeson system is 
produced predominantly in spin states J — (S'-wave) and J — 1 (P-wave). The 
results enable us to study the pion production and the resonance production on 
the level of spin-dependent amplitudes rather than spin-averaged cross-sections. In 
particular, a model independent separation of 5'-wave and P-wave amplitudes is 
possible only in measurement on transversely polarized targets. 

The high statistics measurement of n~p — > TT~iT^n at 17.2 GeV/c at CERN-PS 
on unpolarized target^ was later repeated with a transversely polarized target^ ^ 
at the same energy. The resulting model independent amplitude analysis^'^'^ pro- 
vided the first evidence for significant contributions from helicity-nonflip amplitudes 
with Ai exchange quantum numbers {I*^ = 1~, J^*^ = 1+"'") which had long been 
assumed absent. Moreover, the 5- wave partial wave intensity Is showed a clear 
bump in the 750-800 MeV mass region® in one of the two solutions^'^'^ which indi- 
cated the possibility of a new scalar state with a mass near p^. 

Additional information was provided by the first measurement of tt'^u — > 
TT+TT^p reaction on polarized deuteron target at 5.98 and 11.85 GeV/c also done 
at CERN-PS. ^'^^'^^ Amplitude analysis at larger momentum transfer confirmed the 
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evidence for large ^i-exchange contributions and found resonant-like structures in 
the 5"- wave partial-wave intensities.^^ This data also found important t-dependent 
structures in the moduli of P-wave amplitudes within the p° mass region. -"^^ 

In a recent paper^* we focused on the evidence for the scalar state 1 = 
0++(750) coming from all these measurements on polarized targets. We found 
that all four solutions for S-wave partial wave intensity at 5.98 and 11.85 GeV/c 
suggest resonance structure around 750 MeV. Its tt+tt" decay width depends on 
the solution and was estimated to be in the range of 100-250 MeV. The evidence 
for this state is strengthened by the fact that the 5'-wave amplitudes S and S are 
nearly phase degenerate with the dominant resonating P-wave amplitudes L and 
L. It was suggested that the 1 = 0"''"'"(750) state is best understood as the 
lowest- mass scalar gluonium 0~^~^{gg). Our results are in agreement with the S- 
wave partial-wave intensity for 7r"''7r~ — > 7r°7r° estimated-*^^ from the measurement 
of TT+p 7r°7r°A++ at 8.0 GeV/c. 

The data on transversely polarized targets are best analysed in terms of nu- 
cleon transversity amplitudes. There are two /S-wave and six P-wave amplitudes. 
Amplitude analysis expresses analytically^'^^ the eight moduli and six cosines of rel- 
ative phases of nucleon transversity amplitudes in terms of measured spin density 
matrix (SDM) elements. There are two similar solutions. However, in many (m, t) 
bins the solutions are unphysical: typically a cosine has magnitude larger than one 
or the two solutions for moduli are complex conjugate with a small imaginary part. 
In Ref. 12 and 14 we presented the direct analytical results (taking only real part of 
complex solutions). The authors of Ref. 3-6 took these analytical solutions as start- 
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ing values for a minimization program which fitted the measured observables to 
obtain physical values of moduli and cosines. 

The occurrence of unphysical solutions is a major difficulty for all amplitude 
analyses of pion production and the source of uncertainty in the evidence for / = 
0++ (750) resonance. In this paper we investigate this problem using Monte Carlo 
methods. ^^'^^ 

The basic idea of Monte Carlo amplitude analysis is to filter out^^ the unwanted 
unphysical solutions and to determine^^ the range of physical values of moduli and 
cosines of relative phases. To achieve this we randomly varied the input SDM 
elements within their errors, performed the amplitude analysis for each new set of 
the input SDM elements, and retained the resulting moduli and cosines only when 
all of them had physical values in both solutions. The results presented in this 
report are based on 30 000 random variations of input SDM elements in each (m, t) 
bin. The distributions of moduli and cosines define the range of their physical values 
and their average values^^ in each (m, t) bin. 

The results for the P-wave moduli are essentially the same as in the previous 
analysis^"* but the changes for the S-wave are striking. After filtering out the un- 
wanted unphysical solutions, a clear signal for / = 0++(750) state emerges in all 
four solutions for 5'- wave partial wave intensity at all 3 energies. 

The paper is organized as follows. In Section 2 we review the basic formal- 
ism. In Section 3 we describe our Monte Carlo approach to amplitude analysis of 
ttN — > 7r"'"7r~A'" reactions on polarized target. The evidence for the / = 0"'""'"(750) 
state is presented in Section 4. In Section 5 we describe the spin dependence of p° 
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production and discuss its unexpected features. In Section 6 we discuss the assump- 
tions involved in the determination of tttt phase shifts from data on 7r~p ^ iT~n~^n 
on unpolarized target and present tests of the key assumption of absence of Ai- 
exchange in measurements on polarized targets. The paper closes with a summary 
in Section 7. 

2. BASIC FORMALISM 

For invariant masses below 1000 MeV, the dipion system in reactions ttN — > 
Tv^Tv^N is produced predominantly in spin states J — (^-wave) and J — 1 
(P-wave). The experiments on transversely polarized targets then yield 15 spin- 
density-matrix (SDM) elements describing the dipion angular distribution. The 
measured SDM elements are^°'^^ 

Pss + Poo + 2pii, poQ-pii, pi-i (2.1a) 

Repio, Repis, Repos 
pL + P^o + 2p?i, P^o-Pfi, pU (2.1fe) 
Rep\o, Rep\s, RePos 
Im pI_^, ItupIq, Impl^ (2.1c) 

The SDM elements (2.1a) are also measured in experiments on unpolarized targets. 
The observables (2.1b) and (2.1c) are determined by the transverse component 
of target polarization perpendicular and parallel to the scattering plane ttN — > 
(7r"'"7r~)A^, respectively. The SDM elements (2.1) depend on s,t,m where s is the 
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cm. system energy squared, t is the four-momentum transfer squared, and m is the 
'JT^'JT~ invariant mass. There are two linear relations among the matrix elements in 
(2.1): 

Pss + poo + 2pii = 1 (2.2) 

pL + pIo + 2p?i = A 

where A is the polarized target asymmetry. 

The reaction yr+n ti^ti~p is described by pion production amplitudes 
-f^Ap,OA„ is,t,m,$, (f)) where Xp and A„ are helicities of the proton and neutron, 
respectively. The angles 9, describe the direction of tt^ in the tt+tt" rest frame. 
The production amplitudes can be expressed in terms of production amplitudes 
corresponding to definite dipion spin J using an angular expansion 

oo + J 

i?A„OA„ = E E (2J+ l)^iy/A„OA„(«,*,™)^Ao(^)e'^'^ (2.3) 

J=0 A=-J 

where J is the spin and A the helicity of the 7r"'"7r~ dipion system. Our amplitude 
analysis is carried out in the t-channel helicity frame for the 7r"'"7r~ dimeson system. 
The helicities of the initial and final nucleons are always in the s-channel helicity 
frame. 

The "partial-wave" amplitudes H^^^ can be expressed in terms of nucleon 
helicity amplitudes with definite t-channel exchange naturality. In the case when 
the TT+TT" system is produced in the S- and P-wave states we have 

0-^^-0+^^: H^+,o+ = So (2.4a) 

Ho+,0- = ^1 
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0-i : H'o+,o+-Lo (2.46) 

Ho+,0- = 
1 _ iVo ± £/o 

In (2.4), 0~ stands for pion, |^ for nucleon, 0+ for J = dipion state (S-wave), 
and 1~ for J = 1 dipion state (P-wave). At large s, the amplitudes A^o and Ni 
are both, dominated by natural A2 exchange. The amplitudes S^^ Lji^ Un,n =0,1 
are dominated by unnatural exchanges: Ai exchange for n = and tt exchange for 
n = 1. The index n = |An — Ap| is nucleon helicity flip. 

The data on transversely polarized targets are best analysed in terms of nucleon 
transversity amplitudes (NTA's).^'^^'^^ In our kinematic region we work with two 
S'-wave and six P-wave NTA's of definite naturality defined as follows ^'■'^^'■'^^ 

S = iSo + iSi)/V2 , 5 = (5o - iSi)/V2 (2.5) 

L = (Lo + iLi)/V2 , L = (Lo - iLi)/V2 

U={Uo + iUi)/V2 ,U={Uo- iUi)/V2 

N = {No- iNi)/V2 ,N = {No + iNi)/V2 

The amplitudes S, L,U, N and S, L,U,N correspond to recoil nucleon transversity 
"down" and "up", respectively. ■^■'^'■^^ The "up" direction is the direction of normal 
to the scattering plane defined according to the Basel convention by p^^ x p^^r where 
and Ptttt are the incident pion and dimeson momenta in the target nucleon rest 
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frame. The 5'-wave amplitudes S, S and P-wave amplitudes L, L have dimeson 
helicity A = 0. The pairs of amplitudes U, U and N, N are combinations of nucleon 
helicity amplitudes with dimeson helicities A = ±1 and have opposite t-channel- 
exchange naturality. 

We can now express the observables in terms of amplitudes. In our normaliza- 
tion, the integrated cross section E = (Pa/dmdt is given by 

E = J] \Sn\^ + |L„|2 + |[/J2 ^ |^^|2 (2.6) 
n=0,l 

= |5p + \S\'^ + |L|2 + |Lp + + |C7|2 + |7Vp + |]V|2 

The cross section has not been measured in the experiments on polarized targets. 
Consequently, we will work with normalized amplitudes corresponding to 

S=-^^l (2.7) 

Using (2.6), the relations for SDM elements in terms of normalized helicity ampli- 
tudes read as follows. ^'^^ 

Unpolarized SDM elements 



Pss + Poo + 2pii = J2 + l^nP + |t/nP + \Nn\'^ (2.8a) 

n=0 

P00-Pll= l^np- ^(|iVnP + |C^nP) 

n=0,l 

Pl-1= E l(\Nn\' - \Un\') 



2 

n=0,l 



n=0,l 

V2Repu = E MUnS*J 

n=0,l 
n=0,l 
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Polarized SDM elements 



Pis + Poo + 2p?i = 2Im{S^Sl + LoLl + U^Ul + N^N*) 
Poo - Pii = Im{2LoLl - iVoiV* - C/qC/*) 



(2.86) 



pI_^ = Im{NoN* - UoU^) 
V2ReplQ = Im{UoLl -UiL*) 
V2Repl^ = Im{UoSl - UiS*) 
Replg = Im{LoSl - LiSq) 



-Impf_^ = Im{NoU^ + N^U^) (2.8c) 
V2Implo = Im{NoLl + NiL*^) 
V2Impf, = Im{NoSl + N^S^) 

Only the polarization dependent SDM elements measure the nucleon helicity flip- 
nonflip interference. The observables (2.8b) and (2.8c) measure the interference 
between the amplitudes of the same and opposite naturalities, respectively. 

To express the observables in terms of normalized nucleon transversity ampli- 
tudes (2.5), we first introduce partial wave cross-sections cr{A) and partial- wave 
polarizations t{A) defined for amplitudes A = S,L,U,N as 

a{A) ^ \Ao\^ + \Ai\^ ^ \A\^ + \A\^ (2.9) 

t{A) = 2eIm{AoAl) = \A\^ - \A\^ 
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where e = +1 for A = S, L,U and e = — 1 for ^ = N. In our normalizatfon the 
reaction cross-section is 

E^a{S) + a{L) + a{U) + a{N) ^1 (2.10) 

The relations for SDM elements (2.8a) and (2.8b) in terms of normalized nucleon 
transversity amplitudes (2.5) and quantities (2.9) read 

Pss + Poo + 2pii = a{S) + a{L) + a{U) + a{N) (2.11a) 

Poo-pii^<y{L)-h^a{U) + a{N)] 

Pi-i = -\[cj{U)-a{N)] 

Pis + pIo + 2p?i = t{S) + t(L) + t{U) - t{N) (2.116) 

plo-pli=r{L)-\[T{U)-T{N)] 

p\-i = -\[r{U) + T{N)] 

V2Repio = Re{UL* + UL*) (2.12a) 

V2Repu = Re{US* + UL*) 

Repos = Re(LS* + LS*) 

V2ReplQ = Re{UL* -UL*) (2.126) 

V2Repl^ = Re{US* - C7 5*) 

Rep^^ = Re{LS* -LS*) 

The relations (2.11) and (2.12) suggest to introduce new observables which are the 
sum and difference of the SDM elements (2.8a) and (2.8b). Using the notation of 
(2.2), the first group of new observables reads 

ai = ^[1 + ^] = \S\^ + |L|2 + + \N\^ (2.13a) 



2' 
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02 = [(poo - Pii) + (P^o - Pii)] = 2|L|2 - _ lATp 
a3=\pi-i + pl_^] = \W-\U\^ 

04 = ^[-Repio + Replol = 003(7^^^) (2.136) 

v2 

05 = ;^[^epis + i?epfj = \U\\S\ 003(7^^) 

06 = ^[Repos + Repos] = \L\\S\ 003(7^^) 

Similar equations relate the diffeence of SDM elements to amplitudes of opposite 
transversity. The second group of observables is defined as 



ai = 


^-[1-A] = \S\' + \L 




(2.14a) 


02 = 


[(poo - Pii) - (Poo - 


p^i)] = 2|L|2 - |C7|2 - |Ar|2 




03 = 


[pi-i-pLi] = |iVp 


-\W 




0,4 — 


■^[Repio - Replol -- 


^\U\\L\ cos (7^^) 


(2.146) 


05 = 


■^[Repis - Repu] -- 


= Fl|:S|cos(7^^) 




ae = 


\[RepOs - Repls] = 


\L\\S\cosi^^J 





In the Equation (2.13b) and (2.14b) we have introduced explicitly the cosines of 
relative phases between the nucleon transversity amplitudes. 

The SDM elements (2.8c) form the third group of observables ^'^^ which is not 
used in the present amplitude analysis. 

The first group (2.13) involves four moduli \S\'^, |Lp, |f/p and |iVp, and three 
cosines of relative phases cos(7^^), cos(7gj^) and cos{'fj^jj). The second group (2.14) 
involves the same amplitudes but with opposite nucleon transversity. In Ref. 12 we 
derived analytical solution for these amplitudes in terms of observables. For the 
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first group we obtained a cubic equation for \L\'^ = x 

ax"^ + hx^ + cx + d = i) (2.15) 

with coefficients a,b,c,d expressed in terms of observables ai,i = 1,2, ...,6 (see 
Ref. 1, 12). The remaining moduh and the cosines are given by expressions 

= + a2) -3|L|2 (2.16) 
|t/|' = |Lp-^(a2 + a3) 
\N\^=\L\^-^{a2-as) 

= \iW\ 

cos(75[;) = ^5^: cos(7g^) = 

The solution for the second group (2.14) is similar. 

The physical solutions of cubic equation (2.15) for |Lp must produce physical 
and normalized moduli and physical values for the cosines 

^<\Af <\, A = L,S,U,N (2.17) 

-1 < cos 7fc < 1, A; = Lt/, SU, SL 

There are similar constraints on the solutions for |Lp in the second group. 

The analytical solution of the cubic equation (2.15) is given in Table 1 of 
Ref. 12. One solution of (2.15) is always negative and it is rejected. The other 
two solutions are generally positive and close. However in a number of (m, t) bins 
we get unphysical values for some cosines and in some cases also negative moduli 
of amplitudes. In some (m, t) bins the mean values of input SDM elements yield 
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complex solutions for |Lp or |Lp or both (with positive real parts). To filter out 
the unwanted unphysical solutions we now turn to Monte Carlo amplitude analysis. 

3. MONTE CARLO AMPLITUDE ANALYSIS 

The origin of the presented Monte Carlo method lies in the problem of er- 
ror combination and propagation when the function of input uncertain variables 
is highly nonlinear. This is our case. The analytical solutions for the amplitude 
X = |Lp from the cubic equation (2.15) are highly nonlinear functions of the in- 
put observables - see Table 1 of Ref. 12. Even if the input spin density matrix 
(SDM) elements have Gaussian distributions, the solutions of the cubic equation 
(2.15) and the amplitudes (2.16) are non-Gaussian distributions as the result of the 
nonlinearity. The question arises how to estimate the errors on the amplitudes. 

In general, when the errors on the input variables are small, one can use a 
linear approximation to error propagation (Ref. 17). This method yields symmet- 
ric errors using la errors as input. This approximation was used in our previous 
analyses. Since the errors on polarized SDM elements are not actually small, this 
approximation is not satisfactory. 

In his review paper (Ref. 17) F. James advocates the use of Monte Carlo method 
as perhaps the only way to calculate the errors in case of nonlinear functions which 
produce non-Gaussian distributions. The method has the added advantage that it 
can separate the physical and unphysical solutions, something the linear approx- 
imation and the method cannot do. The Monte Carlo method was first used 
to calculate errors of analytical solutions of a cubic equation by H. Palka^^ in an 
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amplitude analysis of reactions 7r~p — > K'^K~n and 7r~p — > K^K^n at 63 GeV/c. 

In our Monte Carlo amplitude analysis the input SDM elements were randomly 
varied within their experimental errors in each (m, t) bin and amplitudes were calcu- 
lated for each new selection of SDM elements. The input SDM elements were varied 
independently using a multidimensional random number generator SURAND^^ with 
an initial seed number 80629.0. Since any sequence of random numbers generated 
by SURAND is reproducible, our results are also reproducible. 

Each selection of SDM elements yields two solutions for amplitudes of group 1 
(eqs. (2.13)) 

|L|2, |iV|2 (3.1) 

cos(7l[7), cos(75j^), cos(7g^) 
and two solutions for amplitudes of group 2 (eqs. (2.14) 

|t7|^ |iv|2 (3.2) 

cos(7^^), cos(7^^),cos(7^^) 

In each group the solution was classified as physical only when all 4 moduli and 
all 3 cosines of relative phases had physical values. The selection of SDM elements 
was classified as pass only when all solutions for amplitudes (3.1) and (3.2) were 
physical solutions. Otherwise the selection was classified as fail. 

The Monte Carlo amplitude analysis program was run with a total number of 
SDM elements selections Ntot = 10 000, 20 000 and 30 000 in each (m, t) bin. Each 
selection is classified as pass or fail according to the above criteria. The passing 
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rate, or the ratio Npass/Ntoti was nearly constant in each (m, t) bin as Nfot was 
increased from 10 000 to 20 000 and 30 000 selections. Consequently no further 
increases in Ntot were attempted. However, there are considerable variations of 
Npass from bin to bin. In Figure 1 we show the m-dependence of Npass for the runs 
with Ntot — 30 000 for reactions 7r~p-f — > TT~n~^n at 17.2 GeV/c and 7r+n| — tt+tt"]? 
at 5.98 GeV/c. The results at 11.85 GeV/c are similar. For some reason Npass is 
lowest for m ~ 800 MeV at all 3 energies. At each energy there is one bin which 
produced no physical solution (880-900 MeV bin for 17.2 GeV/c, 360-440 MeV bin 
at 5.98 GeV/c, and 820-870 MeV bin at 11.85 GeV/c). At 17.2 GeV we shaU omit 
also the mass bin 640-660 MeV where we found Npass — 3 only. 

The pass and fail SDM elements are retained in each (m, t) bin to determine 
their distribution, range of values and average value. A priori, the ranges and 
average values are expected to be different for pass and fail selections and could 
depend on SDM elements. However, we found that for all unpolarized SDM elements 
in all (m, t) bins the ranges of values for pass and fail Monte Carlo selections are the 
same and coincide with the ranges given by the original errors. Also surprisingly, 
the average values of pass and fail unpolarized SDM elements are equal to the mean 
values of the input SDM elements. For polarized SDM elements we reach the same 
conclusion for fail Monte Carlo selections. For pass Monte Carlo selections, the 
polarized SDM elements have smaller ranges in a few (m, t) bins with lowest Npass 
values and their average values differ up to 6% from the mean values of input SDM 
elements. The situation for p^^ + Pqq + 2p\-^ is typical and is illustrated in the Table 
1 for 7r~p — > n~n^n at 17.2 GeV/c. 
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The values of moduli and cosines obtained for each pass Monte Carlo selection 
were collected to determine their range and average values in each (m, t) bin. To 
study the distribution of values of moduli and cosines for all Npass selections, the 
moduli were binned in bins of size 0.02 and the cosines in bins of size 0.04. An 
example of such distributions is shown in Figure 2 where we present solutions 1 
and 2 for cos(7^^) at 5.98 GeV/c for m =520-600 MeV and -t =0.2-0.4 (GeV/c).^ 
Notice the typical differences between the distributions for Solutions 1 and 2. The 
distribution for solution 1 yields a range of cos from 0.52 to 1.00 with an average 
value 0.94. The distribution for Solution 2 yields a range of cos7g^ from 0.31 to 
1.00 with an average of 0.64. The differences between the two distributions suggest 
that the mimimization approach used in Ref. 3,4,5 to obtain physical solutions 
may not be dependable. 

The resulting ranges of values and average values for all unnormalized moduli 
and for cosines are shown in Figure 3 for 7r~p Tr^Tr+n at 17.2 GeV/c and in 
Figure 4 for Tr+n — > tt~^tt~p at 5.98 GeV/c. The figures show the mass dependence 
of the two solutions for moduli and cosines for —t =0.005-0.2 (GeV/c)^ at 17.2 
GeV/c and -t =0.2-0.4 (GeV/c)^ at 5.98 GeV/c. The unnormalized moduh l^pE 
and |ApE, A — S, L, U.N are calculated using E = d'^a/dmdt from Ref. 2 at 17.2 
GeV/c and from Ref. 10 at 5.98 and 11.85 GeV/c. A comparison with Figures 
1 and 2 of Ref. 14 reveals that the physical solutions for S'-wave moduli and for 
cosines are much smoother and continuous. The jitter which characterized the 
unphysical solutions is removed by filtering out the unphysical solutions. We discuss 
the features of the amplitudes in detail in the next two sections. 
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For each pass Monte Carlo selection we also calculated the normalized partial- 
wave cross-sections aA and partial-wave polarizations ta defined in (2.9) in order 
to obtain their range of values and average values. Using experimental results for 
S = d^a/dmdt, we then calculated partial wave intensities I a = caS where A = 
S, L, U, N. The solutions for amplitudes with opposite transversities are entirely 
independent. We can denote the two solutions for normalized nucleon transversity 
amplitudes as A{i) and A{j) with i = 1,2 and j = 1, 2. Because the moduli squared 
in (2.9) are independent, there is a fourfold ambiguity in the partial-wave intensities. 
Using the indices i and j to label the four solutions, we get 

lA{t,j) = [\A{t)\' + \A{j)f]j: (3.3) 

where A = S,L, U, N. The partial-wave intensities (3.3) are defined in the physical 
region and depend on the dipion mass m and momentum transfer t. 

The Monte Carlo amplitude analysis presented in this paper involves one im- 
portant simplification, namely the experimental data points are represented by a 
rectangular distribution of width 2a. The 2a width is motivated by the requirement 
to study the propagation of la errors on experimental data through the highly non- 
linear expressions for the amplitudes in order to define the errors on amplitudes. 
The uniformity of the experimental data distributions is due to the use of the uni- 
form random number generator. 

The actual experimental data are presumed to have Gaussian distributions. To 
examine the effect of Gaussian distribution of data on the results of amplitude anal- 
ysis, we analysed the data at 17.2 GeV/c and 5.98 GeV/c, using an IBM Gaussian 
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random number generator SNRAND.-^^ 

In this calculation the "fail" distributions of SDM elements are approximately 
Gaussian with a broad range of half-size of 4.-4.5(7. The "pass" distributions of 
SDM elements (giving the physical solutions for amplitudes) have a narrower range 
of approximate half-size 2cr. Most of the points lie within la half-interval. Between 
1(7 and 2(7 the "pass" distributions have only low tails. This means that most 
physical solutions lie within la vicinity of central points and none beyond 2a. 

As expected, the distributions for amplitudes (moduli and cosines) are non- 
Gaussian. Importantly, the average values of distributions are essentially identical 
to the average values obtained previously. The agreement is such that the black 
points in the figures do not change. The range of distributions is now broader. The 
minimum and maximum values change typically by 0-20%. However this increase is 
due to only low tails of distributions. The broadening of range of values now makes 
more difficult the estimation of errors. The usual calculation of standard deviation 
does not make sense in the case of such highly non-Gaussian distribution as is e.g. 
Solution 1 for COS75L (Fig. 2). 

Because the average values of amplitudes remain virtually the same and the 
distributions for amplitudes acquire only low tails, there is no significant change 
from the results presented in this paper, and the simplification used in our Monte 
Carlo amplitude analysis is an acceptable approximation. 

4. EVIDENCE FOR A SCALAR STATE 7 = 0++(750) 

After filtering out the unphysical solutions, the evidence for the scalar state / = 
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0"'"'"(750) becomes much more convincing than were the previous indications.^^ 
The evidence is based on three observations: 

(a) Resonance structure of the 5'- wave partial wave intensity Is- 

(b) The resonant structure of unnormahzed ^-wave moduh. 

(c) The constant relative phase between the S'-wave amplitudes S, S and the dom- 
inant P-wave amplitudes L, L. 

We first look at the S-wave moduli at 17.2 GeV/c at lower momentum transfers 
which are summarized in Fig. 5. We see that the amplitude \S\'^Tj has a clear 
resonant structure in both solutions. This is a change from previous results in 
Ref. 14 where only the solution 1 was clearly resonating. In both solutions the 
maximum is at 750-770 MeV with a width at half-height estimated at 175-200 
MeV. The moduli \S\'^'Z of opposite transversity have smaller magnitudes and do 
not show a resonant structure. At first sight this might be surprising, but we see the 
same effect also in p° production amplitudes. For instance, while the amplitudes 
|LpS, iC^pS and lA'^pS show resonant structures, the amplitudes |f/pE and 

|A^|^S are smaller and show less resonant structure. We can ascribe this effect to the 
spin dependent dynamics of resonance production. The difference of magnitudes of 
\S\'^Tj and \S\'^Tj is due to the presence of Ai-exchange i.e. nonvanishing amplitude 
-S'o. 

The usual signature of a production resonance is a peak or bump in the pro- 
duction cross-section. The S-wave partial-wave intensity Is = {\S\^ + \S\'^)T, at 17.2 
GeV/c is shown in Fig. 6. We find a clear resonant signal at or near 750 MeV in 
two solutions /^(l, 1) and /s(2, 1). The other two solutions /s(l, 2) and 15(2, 2) do 
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not show a clear resonant structure. In these solutions Is is rising up to 770 MeV, 
then it drops but remains high. The cause of this behaviour is the nonresonating 
amplitude \S\'^T, which is large above 800 MeV. We are thus faced with an apparent 
ambiguity in the 5'-wave intensity observed also in the earlier studies. ^'^'^ However 
more fundamental evidence for a"(750) production at 17.2 GeV/c at lower t is the 
fact that both solutions for the amplitude resonate. We see the ambiguity in 

Is only because the amplitudes |<S'pS and |<S'pS happen to have magnitudes which 
are not very different. We do not see any ambiguity e.g. in the intensity II because 
the amplitude l-C/pS is much larger than jl/pS which also remains high above 800 
MeV. 

As seen in Figure 3, the amplitudes l-C/pS and \L\'^T, are dominant p° production 
amplitudes. The relative phase between the (5- wave amplitudes S and S and the 
P-wave amplitudes L and L is another important piece of evidence for existence of 
the 1 = O"*""*" (750) resonance. In Solution 1 both relative phases and 7^^ are 
consistent with zero and we have a phase degeneracy of amplitudes S and L, and 
S and L. In Solution 2 the relative phases are not zero but are small and constant 
over the considered mass region. Since the amplitude L is resonating, the phase 
degeneracy with amplitude S suggests that S also resonates with a resonance mass 
near po- 

The S'-wave intensity at larger momentum transfers —t = 0.2 — 0.4 (GeV/c)^ at 
5.98 and 11.85 GeV/c is shown in Fig. 7 and 8, respectively. There is no ambiguity 
at larger momentum transfers. At both energies we find a clear resonant signal at or 
near 750 MeV in all solutions. At 5.98 GeV/c, the width at half-height is estimated 
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to be 270 MeV. We note that at 11.85 GeV/c, the solutions Is{2, 1) and /s(2,2) 
are consistent with a narrow width of 150 MeV. 

The resonant structure of moduh at larger momentum transfers —t =0.2-0.4 
(GeV/c)^ at 5.98 GeV/c and 11.85 GeV/c is less clear. At both energies it is the 
Solution 2 which has a stronger indication for a resonance in both moduli 
and j/SpS. Because of lower statistics, the errors on input SDM elements are larger 
which results in larger ranges of values for cosines of relative phases. Nevertheless, 
we also observe that the relative phases and 7^^ are near zero in Solution 1, 
and are only slowly varying over the critical mass region in Solution 2. The 
unusually larger variation of cos(7^^), cos(7^^), cos(7^^) and cos(7^j^) at m =980 
MeV could be due to the presence of scalar resonance /o(975). 

We have also fitted the S'-wave partial-wave intensities Is to a Breit-Wigner 
form using the CERN optimization program FUMILI.^-*^ The parametrization used 
had the general form 

lAim) ^ Na\BWr\'' (4.1) 

where A — S, L,U, N and 

BWn = (^)V2J+T , "^f. ^ (4.2) 

In (4.1) A^^ is a normalization constant which includes square of elasticity x and 
isospin factor. In (4.2) q is the 7r~ momentum in the tt^tt"'" rest frame 



q=\q\ = V0.25m2 - (43) 

where /x is the pion mass. In (4.2), J, niR and F are spin, mass and width of the 
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resonance. The mass dependent width is 

r-r,(^)--^^f^ (4.4) 

Dj{qr) 

where qn = q{m = mn) and Dj are centrifugal barrier functions of Blatt-Weiskopf^^ 

Doiqr) = 1. (4.5) 
Di{qr) = 1. + {qrf 

In (4.5) r is the interaction radius. The factor mj ^Jq in (4.2) arises from the Chew- 
Low formula. 

For S'-wave ^ = 5', J = 0, i? = cr and F is independent of r. The results of 
the fit to the two resonating solutions for Is at 17.2 GeV/c are presented in Table 
2. We encountered difficulties in fitting the form (4.1) to the solutions /s(l, 2) and 
75(2, 2). As seen in Fig. 6, the maximum values of the fits (4.1) are systematically 
well below the maximum values of the experimental points for Ig- We conclude 
that the parametrization (4.1) does not represent the experimental data very well. 
The likely reason for this is the sudden drop of values of Is at 790 and 810 MeV 
which can be traced back to the dip in jS'pE at 790 MeV (see Fig. 5). 

The results of the fits to the S'-wave intensity Is at 5.98 GeV/c are presented 
in Table 3 and shown in Figure 7. Again, the maximum values of the Breit-Wigner 
fits (4.1) are systematically below the maximum experimental points of Is- Because 
of large errors on Is at this energy, the resonance parameters in Table 3 have also 
larger errors. 
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5. SPIN DEPENDENCE OF pO PRODUCTION 

Resolution of the p° peak seen in the reaction cross-section dPa/dmdt into its 
individual spin components (8 moduli of nucleon transversity amplitudes) is shown 
in Figures 3 and 4 at 17.2 and 5.98 GeV/c, respectively. We stress that this model 
independent resolution is possible only in measurements on polarized targets. We 
discussed the S-wave contributions in the preceeding Section. In this Section we 
focus on the P-wave contributions. We observe several important features in the 
spin dependence of production on the level of amplitudes. 

(a) The mass spectrum on the level of amplitudes depends on dimeson helicity A 
and on the spin of the recoil nucleon. At low momentum transfer (Fig. 3) the 
unnatural exchange amplitudes and |t/pS with recoil nucleon spin up are 
larger than the amplitudes |ivpS and |C/pS with recoil nucleon spin down. The 
opposite is true for natural exchange amplitudes |A^pE and |A^pE. These findings 
are similar to observations reported in Fig. 10 of Ref. 4. At larger momentum 
transfers (Fig. 4) we also find \L\^T, and |A'"pS larger than |LpS and |A'"pS, 
respectively, but |t/pE ^ These large differences between the unnatural 
amplitudes with recoil nucleon spin up and down are due to the contributions from 
large and nontrivial Ai-exchange amplitudes Lq and Uq. The same interpretation 
was given in Ref. 4 and 5. We note that the t-dependence of amplitudes in the p° 
mass region also shows evidence for large Ai exchange amplitudes (see Fig. 7 and 
15 of Ref. 4, and Fig. 1 and 3 of Ref. 12). 

(b) We find unexpected suppression of production in several amplitudes. At low 
momentum transfer (Fig. 3) the moduli |A^pS and |C/pS are small and flat and do 
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not show the expected p° peak. Also, while the amplitude jl/pS is still large, it does 
not show a clear peak but rather a broad structure. At the larger momentum 
transfers (Fig. 4) we find p° suppression in the amplitudes |t/pS and |t/pS which 
are small relative to the other P-wave amplitudes. 

(c) The mass spectra of the dominant amplitudes and show unexpected 

narrow structures within the p° mass region which are not seen in the spin-averaged 
cross-section II. The P-wave intensities II, In and lu are shown in Fig. 9 and 10 
for 5.98 and 17.2 GeV/c, respectively. We notice the narrow range of values in these 
partial- wave intensities compared to the range of the moduli. This indicates that 
the values of moduli with recoil nucleon spin up and down are highly correlated for 
each P-wave amplitude. 

First we look at the amplitudes at 5.98 GeV/c (Fig. 9). The mass spectrum of 
the recoil nucleon spin up amplitude shows a clear dip at 757 MeV and a peak 

at 807 MeV. The opposite behaviour is seen in the nucleon spin down amplitude 
|LpE which peaks at 757 MeV and has a dip at 807 MeV. These spin correlated 
structures within po mass region do not appear in the partial-wave intensity II 
which shows a structureless p° line shape. 

Next we note a similar situation in the natural exchange amplitudes but at a 
different mass. The spectrum of jATpS shows a maximum at 782 MeV which is 
associated with a pronounced dip in |A^pE at the same mass. The partial wave 
intensity In again shows a structureless line shape expected from a p° 

One should note however that the apparently most significant (2 — Scr) narrow 
structure at 807 MeV is due to a single deviating date point in the p\_i — (see 
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Fig. 8 of Ref. 10). This deviation is not observed in the p^^ + pgo + '^Pii^ ^^^^ 
involving the difference between |Lp and |Lp. This shows that the structures at 
5.98 GeV/c can be treated as indications only. 

Partial- wave intensity II at 17.2 GeV/c at lower momentum transfer (Fig. 10) 
shows a structureless line shape which peaks at 790 MeV. The expectation that the 
same structureless line shape is reproduced on the level of amplitudes is not met 
again. Instead we find spin correlated structures in the amplitude and 
The amplitude [L^E peaks at 790 MeV and has a lower value at 770 MeV. The 
opposite is seen in the amplitude l-C^pS which has a lower value at 790 MeV and 
peaks at 770 MeV. 

The partial- wave intensities presented in Fig. 9 and 10 are solution 
/a(1, 1), A = L^U, N. The other solutions are similar. 

At this point one may ask how significant is the evidence for the structures in 
the p line shape on the level of amplitudes. We note that the structures occur in 
(m, t) bins with relatively high statistics which provides a measure of confidence in 
the results. Also, the t-evolution of mass dependence of moduli of normalized am- 
plitudes shows t-dependent structures within p° mass region (see Fig. 1 of Ref. 13). 
Another possibility is to fit Breit-Wigner form (4.1) to the spin-averaged partial 
wave intensity II and superimpose this Breit-Wigner fit (suitably normalized) over 
the mass spectra in moduli and The numerical results of the fits 

to II are given in Tables 4 and 5. In Fig. 11 we show the results for reaction 
TT+n TT+7T-p at 5.98 GeV/c and -t = 0.2 - 0.4 (GeV/c).^ The Breit-Wigner 
fit to II peaks at 780 MeV. Superimposed over the moduli and |-^pE, the 
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Breit-Wigner fit shows that the narrow structures in their mass spectra are sta- 
tistically significant. At 757 MeV, the structures are la effect (dip in |LpE and 
peak in |LpE). At 807 MeV, the structures are 2-3cr effect (peak in \L\'^T, and 
dip in |LpE). In Fig. 10 we show the results for 7r~p — > tt^tt+p at 17.2 GeV/c 
and —t — 0.005 — 0.2 (GeV/c).^ While the Breit-Wigner form fits the intensity II 
very well we see that the mass spectrum of amplitude jl/pE is narrower than the 
Breit-Wigner fit and the mass spectrum of amplitude is broader than the 

Breit-Wigner fit to II- At present, the evidence for structures in the p-line shape 
of amplitudes at larger t is still preliminary and insufficient for theoretical analyses. 
New experiments on polarized targets with significantly higher statistics than the 
60,000 events (at 5.98 GeV/c) in the Saclay measurement -"^^ are required to confirm 
the existence of such structures and to investigate the t-dependence of p line shape 
on the level of amplitudes. 

We now briefiy discuss the physical significance of the observed spin effects in 
production. The presence of Ai-exchange (feature (a)) will be discussed in the 
next Section. Here we will focus on the features (b) and (c). 

Following the discovery of resonances in hadron interactions, it has always 
been believed that the production and decay of resonances were separate events. 
For instance, the reaction ir'^n — > 7r"''7r~p was thought of as a two step process in 
the p^ mass region: p^ production Tr+n —>■ p^p followed by p° decay p^ —>■ tt+tt". 
In this picture the ^-matrix elements factorize into production and decay matrix 
elements 

T{7i+n TT+TT-p) = T{n+n p^p)(t){p^)T{p^ Tr+Tr") (5.1) 
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where 0(p°) is a p° propagator leading to the Breit-Wigner description of dipion 
mass dependence of modulus of each production amplitude \T{TT~^n —>■ iT~^n~p)\'^ as 
well as the production cross-section^^ E = a/dmdt. 

It is expected from (5.1) that the same p° resonance line shape seen in the spin- 
averaged cross section d^ a/dmdt will appear also in the modulus of every P-wave 
spin-dependent production amplitude. The suppression of p° production observed 
in several spin amplitudes and described above in part (b) seems to invalidate the 
factorization hypothesis and the underlying simple picture of resonance production. 
However, there is no theoretical explanation for the suppression of p° production 
in these spin amplitudes (Fig. 3 and 4). Evidently, the spin amplitudes showing po 
suppression cannot be fitted with Breit-Wigner form. 

The factorization hypothesis (5.1) also implies that the line shape of the mass 
spectrum of the decaying resonance does not depend on the nucleon spin and on the 
momentum transfer t. However, the line shapes of the dominant amplitudes jl^pS 
and at larger t show unexpected structures within the mass region which 

correlate the mass spectra corresponding to opposite nucleon spins. Comparison 
of the structures in Fig. 9 and 10 indicates a change as we go from low to larger 
momentum transfer. This suggests that the narrow dips and peaks observed in 
these moduli within the p° mass region evolve with t (see also Ref. 13). This t- 
dependence of the resonance line shape in the amplitudes is also at variance with 
the factorization hypothesis. 

The suppression of p° production and the narrow structures observed in domi- 
nant amplitudes in mass region represent a new information on pion production 
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and raise questions about the nature of hadron resonances. A resonance which is 
a pole in the amphtude would lead to simple resonance peaks in all moduli with- 
out structures within their widths. We also note that the usual models^^ of meson 
resonances as qq bound states do not predict their hadronic widths and line shape. 

6. TESTS OF ASSUMPTIONS USED IN DETERMINATIONS OF tttt 
PHASE SHIFTS. 

We have presented a model independent and solution independent evidence for 
the scalar state / = 0"''"'"(750) in measurements of ttN-^ — > 'k'^tt~N on polarized 
target. The question arises how to understand the absence of such a state in the 
accepted solution for S-wave phase shift 5q in tttt scattering.^'^^"^^ 

Of course, there are no actual measurements of pion-pion scattering and there 
is no partial-wave analysis of tttt — > tttt in the usual sense. The tttt phase shifts 
are determined indirectly from measurements of 7r~p — > 7r~7r"''n on unpolarized 
target using extrapolations into unphysical region of momentum transfer t and 
several necessary enabling assumptions. Some of these crucial assumptions lead to 
predictions for polarized spin-density-matrix (SDM) elements and are thus directly 
testable in measurements on polarized targets. As we shall see in detail below, these 
assumptions are invalidated in a major way by the data on polarized targets. We 
must use the results of measurements on polarized targets to judge the validity of tttt 
phase shift determinations, and not vice versa. We are thus led to the conclusion, 
that the indirect and model dependent determinations of tttt phase shifts cannot be 
correct. This explains the absence of / = 0"'"+(750) resonance in Sq phase shift 
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from these analyses. 

To unveil the assumptions used in determinations of tttt phase shifts we will 
trace step by step the procedure used by Estabrooks and Martin. ^^'^^ Their ap- 
proach has the advantage of being the most transparent and of using the least 
number of assumptions. Moreover, their assumptions are common to all other de- 
terminations of TTTT phase shifts from unpolarized ttN — > tt+tt^A?" data.^'^^'^^ We 
shall restrict our discussion to dipion masses below 1000 MeV where S- and P-wave 
dominate. 

The starting point are dimeson helicity A = pion exchange amplitudes Si and 
Li in the t-channel. It is assumed that the t- and m-dependence in these amplitudes 
factorizes^'^'^^ 



Si{m,t) = -^Fo{t)—Mm) (6.1) 
Li(m,t) = ^^Fi(t)^v^/i(m) 

where t is the momentum transfer at the nucleon vertex, m and q are dipion mass 
and 7r~ momentum in the tt^tt"'' cm. frame. The form factors Fj{t) describe 
the t-dependence and the functions fj{m), J = 0, 1 describe the mass dependence. 
Further, the functions fj{m) are assumed to be partial wave amplitudes in tt^tt"'' — > 
TT^TT"*" reaction at cm. energy m 

fo - Iff + Iff (6.2) 
/i = ff 

The partial wave amplitudes /J with definite isospin / are defined so that in the 
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TTTT elastic region 

/i = sin55e^^^ (6.3) 

The TTTT phase shifts 5j can be obtained by extrapolating the production amplitudes 
Si and Li from the physical region t < to t — /j,^. One cannot determine both 
1 = and 1 = 2 S-wave phase shifts and so values for /q obtained in analyses of 
TT+p — > TT+TT+n data were used.^'^'^^ 

There is no theoretical proof of factorization (6.1) and identification (6.2) of 
functions fj with tttt partial wave amplitudes. Strictly speaking, the relations (6.1)- 
(6.3) are definitions of tttt phase shifts. It is by no means obvious that these phase 
shifts would coincide with tttt phase shifts determined directly from real pion-pion 
scattering. Only such comparison would test the assumption (6.2) which is not 
testable in measurements of ttN — > Tr+Tr^AT. 

Essential element in the determination of tttt phase shifts is the knowledge of 
the production amplitudes and Li. But an inspection of equations (2.8a) reveals 
that the amplitudes Si and Li cannot be calculated from the data on unpolarized 
target without additional assumptions. There are simply more amplitudes than 
data. To proceed further all determinations of tttt phase shifts must assume that 
all ^i-exchange amplitudes vanish, i.e. 

So = Lo = Uo = (6.4) 

With the assumptions (6.4) and notation (2.9) the equations (2.8a) now read 

Cl = pss+P00 + 2pii = \Sif + \Lif + + Oriv = 1 (6.5) 
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C2 = Poo - Pll = \Li\^ - ^\Uif - ^CTJV 
C3 = Pl-1 = 2'^-'v ~ 2' ' 

C4 = V2Repio = Re{UiLl) = \Ui\\Li\ cos{xlu) 
C5 = y/2Repu = Re{U^Sl) = \U^\\Si\ cos{xsu) 
cq = Repos = Re{LiSl) = |Li | l^i | cos{xsl) 
where the relative phases satisfy identity 

^LU - ^SU + ^SL = (^Li - <^[7i) - (<^Si - <^C/i) + (<^Si " = ^ (6-6) 

The equations (6.5) are formally similar to the set (2.13) and can be similarly solved. 
We obtain 



\Si\' = 


1 + 2C2 




l^il'- 


aN = 








l^lllt^ll 


cos Xsu = 




\Si\\U,\ 



|2 



(6.7) 



- ^6 

For |Lip = a; we have a cubic equation 

„3 , 1^2 



where 



ax'^ + bx'^ + cx + d^O (6.8) 

a = 3 
6 =(1 + 5c2 + 3c3) 

c =(1 + 2c2)(c2 + cs) - 3cl + 4 + 4 

d =(1 + 2C2)C4 + (C2 + cs)4 - '^C^C^Cq 
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Using the data on unpolarized target at 17.2 GeV/c [Ref. 2], Estabrooks and 
Martin found^^'^^ in all (m, t) bins 3 real solutions. One of the solutions is always 
negative and it is rejected. The other two solutions for yield two solutions 

for l^ip. Thus there are two solutions for S- and P-wave phase shifts below 1000 
MeV. Solution 1 for Sq is nonresonating while the solution 2 resonates at the mass 
around 770 MeV with a width about 150 MeV. The resonating solution was rejected 
because it disagreed with a tt^tt^ mass spectrum from a low statistics experiments^ 
on 7r~p — > 7r°7r°n at 8 GeV/c. 

The absence of Ai exchange amplitudes is absolutely crucial for the determina- 
tion of rni phase shifts from unpolarized target data on n~p n~n~^n. Fortunately, 
it is an assumption that can be tested directly in experiments on polarized targets. 
Using the definitions (2.4) of nucleon transversity amplitudes, the assumptions (6.4) 
imply that the moduli of unnatural exchange amplitudes with recoil nucleon spin 
up and down must be equal, i.e. 

\S\^=\S\^, \L\'^ = \Lf, \U\'^ = \U\^ (6.9) 

These predictions are in sharp disagreement with the results of model independent 
amplitude analysis of the data on polarized target at 17.2 GeV/c (Fig. 3) which 
finds 

\S\^ > |L|2 > |L|^ \Uf > \U\^ (6.10) 

Particularly important is the large difference between |Lp and |Lp (see also Fig. 10) 
which can be accounted for only by a large and nontrivial yli-exchange contribution 
from the amplitude Lq. The assumptions (6.4) are thus badly violated. This implies 
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that the determination of amphtudes l-C/ip and \Si\'^ through equations (6.7) and 
(6.8) cannot be correct. This in turn means through (6.1) that the tttt phase shifts 
cannot be correctly determined from the data on unpolarized target. 

Using the assumptions (6.4) in equations (2.8b) leads to the following predic- 
tions for polarized SDM elements: 

pL + pIo + 2p^i = -2(p^o - Pli) = +2pLi (6.11) 

Rep\o = Repl = Repls = (6-12) 

The data for polarized SDM elements clearly rule out these predictions as is shown 
in Fig. 12 and 13 for tt~p —>■ n~TT~^n reaction at 17.2 GeV/c. We find that p^g + 
Poo + 2pii and — 2(pqq — p^^) have large magnitudes but opposite signs while 2p^_^ 
has a small magnitude. The interference terms RepiQ, Repfg and RepQg are all 
dissimilar and have large non zero values. On the basis of this evidence we again 
must conclude that the determinations of tttt phase shifts from unpolarized data on 
ttN — > TT~^TT~N are questionable. 

If the pion exchange amplitudes and Li cannot be reliably determined from 
the data on unpolarized targets, the question arises whether these amplitudes can be 
determined from the data on polarized targets. Unfortunately the data on polarized 
target do not allow the separation of tt- and ^i-exchange amplitudes. In fact it 
has been shown recently^^ that the pion exchange amplitudes Si and Li can be 
expressed in terms of data on polarized target plus the ^i-exchange amplitudes 
^0 and Lq. These relations show explicitly that the determination of tttt phase 
shifts from data on polarized target depends on the model used for Ai-exchange 
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amplitudes and Lq. We conclude that in the absence of a reliable model for Ai- 
exchange or real pion-pion scattering data the tttt phase shifts remain undetermined. 

It is of interest to compare our results for S-wave intensity with previous tttt 
phase shift analyses to assess the effect of exchange and the contribution of 
isospin 1 — 2 S'-wave amplitudes. 

In Fig. 14 and 15 we compare S'-wave cross-sections normalized to one at max- 
imum value. The data points correspond to /s(l, 1) and /s(2, 2) at 17.2 GeV/c in 
the physical region of t. The curves in Fig. 14 (taken from Ref. 15) are predictions 
of various determinations of tttt phase shifts. The curves A and D are the two 
solutions of Estabrooks and Martin (Ref. 27). The curves B and C are predictions 
from Ref. 32 and 33, respectively. The predictions A, B and C show broad struc- 
tures in Is around 750 MeV where the data on polarized target require a narrower 
structure. The broad structure of predictions A^ B and C does not agree well with 
the current algebra prediction above threshold (curve Ref. 34). 

The solid curves in Fig. 15 are prediction from a recent fit to 5q phase shift 
data by Zou and Bugg.^^ Results from another recent fit by Au, Morgan and 
Pennington^^ are very similar. The predictions are broad and peak at 860 MeV 
while our data are much narrower and peak at 750-770 MeV. The predictions from 
phase shift analyses^^'^^ were calculated using Is = |/oP where /o = |/o^° with 
/o"° given by (6.3). 

All curves in Fig. 14 and the solid curves in Fig. 15 correspond to / = 5'=wave 
intensity in tt+tt" — > tt'^tt~ reaction. The data points in Fig. 14 and 15 correspond 
to /S-wave intensity in 7r~p — > 7r~7r "•"??. at 17.2 GeV/c in the physical region of t and 
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include contribution from isospin 1 = 2 S-wave amplitudes. The question arises to 
what extend the 1 = 2 S'-wave amplitudes affect the resonance shape of 5(750) state 
seen in the S'-wave intensity data 7s(l, 1). We cannot answer this question using 
the data on polarized target since these data do not allow isolation of the two 1 = 2 
S'-wave amplitudes. However, we can use (6.2) and the data for / = 2 phase shift 
5q to calculate the S-wave intensity Is = |/oP in tt+tt" scattering with the 1 = 2 
contribution. 

The experimental information about 1 = 2 S-wave phase shifts comes from the 
analysis of Tr+p Tr+Tr+n reaction. We used the results from the 1977 analysis 
of Hoogland et al. in Ref. 37. This analysis presents two very similar results for 
5q based on two different methods (method A and method B). The dashed line in 
Fig. 15 shows the S'-wave intensity of Zou and Bugg with the 1 = 2 correction using 
5q from method A in Ref. 37. The method B yields essentially identical curve. We 
see that the 1 = 2 correction maintains the overall shape of S-wave intensity in 
7r"''7r~ — > TT+TT" and does not secure agreement with the data from 7r~p — > Tr^Tr+n 
analysis. We can conclude that the principal cause of the difference between the tttt 
phase shift analyses and our amplitude analysis of iT~p — > iT~iT~^n data on polarized 
target is the existence of large Ai exchange amplitudes in this reaction. 
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7. SUMMARY 

The measurements of reactions 7r~p-i — > tt^tt+tt- at 17.2 GeV/c and Tr+n-i- — > 
TT^TT~p at 5.98 and 11.85 GeV/c on polarized targets enable a model-independent 
amplitude analysis for dipion masses below 1000 MeV where S-wave and P-wave 
contributions dominate the pion production process. The amplitude analysis yields 
two similar solutions for 8 moduli and 6 cosines of relative phases. In most (m, t) 
bins the analytical solution produces unphysical values of some cosines or moduli, 
violating the conditions (2.17). These conditions are equivalent to inequalities con- 
straints on SDM elements. -"^^'^^ To avoid unphysical solutions the constraints should 
be imposed during the optimization of maximum likelihood function in data tapes 
analysis in future experiments. This will require the use of methods for constrained 
optimization. i6,39,4o,4i 

The previous amplitude analyses^'^ found evidence for a new scalar state / = 
0"''"'"(750) which can be interpreted as the lowest mass scalar gluonium 0'^'^{gg).^'^ 
Because of the significance of this finding it is important to see what happens to 
the scalar mass spectra when the contamination by unphysical solutions is removed. 
To filter out the unphysical solutions we used a Monte Carlo approach. The input 
SDM elements were randomly varied within errors 30 000 times in each (m, t) bin 
and amplitude analysis was performed for each selection. Only when the 8 moduli 
and 6 cosines attain physical values in both solutions are their values retained 
and partial-wave intensities and polarizations are calculated. All such physical 
solutions are collected and these collections then define the distribution, range and 
average values for each modulus, cosine of relative phase, partial wave intensity and 
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polarization. 

After filtering out the unwanted unphysical solutions, a clear signal for / = 
O"*""*" (750) resonance emerges in all solutions for partial wave intensity Is at larger 
momentum transfers at 5.98 and 11.85 GeV/c, and in two solutions at lower mo- 
mentum transfers at 17.2 GeV/c. The ambiguity at 17.2 GeV/c is caused by the 
non resonating amplitude \S\'^T, which masks the resonating behaviour of 
The mass spectra in Is peak at or near 750 MeV. The width at half-height is about 
230 MeV and 270 MeV at 17.2 and 5.98 GeV/c, respectively. At 11.85 GeV/c, 
the solutions 15(2, 1) and 15(2, 2) are consistent with a narrow width of 150 MeV. 
At 17.2 GeV/c, the amplitude \S\'^T, clearly resonates in both solutions at 750- 
770 MeV with a width at half- height estimated to be 175 MeV. The production of 
/ = 0"'""'"(750) appears suppressed in the amplitude IS'pS. There is a clear phase 
degeneracy between amplitudes S and L, and S and L, also indicating resonant 
structure of the S-wave near p° mass. 

The results for p° production generally confirm the findings of previous 
analyses. We find significant suppression of p° production in amplitudes 
|A'"pE and |f/pE at all energies. There is an additional p° suppression in \L\'^T, at 
17.2 GeV/c and in at 5.98 GeV/c. Moreover, the line shapes of and 

|LpS show unexpected structures within p° mass region which correlate the mass 
spectra corresponding to opposite nucleon spins. These narrow structures are not 
observed in the spin averaged partial wave intensities which show the expected 
structureless peak. There is no theoretical explanation for these important fea- 
tures of p*^ production which contradict the factorization hypothesis. 
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All determinations of tttt phase shifts from unpolarized data on tt~p — > Tr^Tr+n 
rely on the absence of Ai-exchange amplitudes. This essential assumption leads to 
predictions for polarized SDM elements which are clearly ruled out by the data on 
polarized target. In our Monte Carlo amplitude analysis we find a clear evidence 
for large and nontrivial Ai-exchange contributions from S'-wave and P-wave ampli- 
tudes, in particular at 17.2 GeV/c. On the basis of this evidence we come to the 
conclusion that the usual determinations of nn phase shifts cannot be correct. This 
also explains why the accepted solution for 5q phase shift shows no evidence for the 
(t(750) resonance. 

Production of scalar state 7 = 0"'""'"(750), the suppression of produc- 
tion in certain production amplitudes and the possible narrow structures in mod- 
uli in mass region are unexpected and important findings of measurements of 
ttN — > 'K'^TT~N reactions on polarized targets. These new phenomena should be 
further investigated experimentally in a new generation of dedicated high statistics 
experiments with polarized targets. ^^'^^ 
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INPUT PASS FAIL 

MASS MEAN ERROR MIN MAX AVRG MIN MAX AVRG 

0.610 -0.702 0.096 -0.796 -0.606 -0.653 -0.798 -0.606 -0.703 

0.630 -0.620 0.085 -0.704 -0.535 -0.585 -0.705 -0.535 -0.624 

0.650 -0.762 0.078 -0.823 -0.808 -0.816 -0.840 -0.684 -0.761 

0.670 -0.596 0.067 -0.663 -0.529 -0.600 -0.663 -0.528 -0.594 

0.690 -0.535 0.057 -0.584 -0.479 -0.509 -0.592 -0.479 -0.534 

0.710 -0.592 0.050 -0.642 -0.542 -0.585 -0.642 -0.542 -0.591 

0.730 -0.610 0.046 -0.656 -0.564 -0.595 -0.656 -0.564 -0.610 

0.750 -0.588 0.043 -0.631 -0.547 -0.591 -0.631 -0.546 -0.587 

0.770 -0.489 0.043 -0.532 -0.447 -0.486 -0.532 -0.447 -0.488 

0.790 -0.557 0.043 -0.598 -0.514 -0.543 -0.599 -0.514 -0.555 

0.810 -0.411 0.046 -0.410 -0.366 -0.380 -0.457 -0.365 -0.410 

0.830 -0.436 0.053 -0.489 -0.383 -0.420 -0.489 -0.383 -0.436 

0.850 -0.461 0.057 -0.488 -0.404 -0.426 -0.518 -0.404 -0.460 

0.870 -0.340 0.060 -0.399 -0.280 -0.315 -0.401 -0.280 -0.342 

0.890 -0.372 0.064 _ _ _ .0.436 -0.308 -0.371 



Table 1. The mass dependence of p^^ + Pqq + 2p\^ in 7i~p — > TT~'K'^n at 17.2 GeV/c 
for —t =0.005-0.2 (GeV/c).^ The input mean values and errors are compared with 
the ranges and average values of p^^ -\- Pqq -I- '2p\i for pass and fail Monte Carlo 
selections for Ntot =30 000. 



45 



Is 


ma- 




Ns 


xVd.o.f. 


Solution 


(MeV) 


(MeV) 






(1,1) 


764± 6 


273 ± 24 


1.97 ± 0.07 


0.360 


(2,1) 


761± 12 


290 ± 54 


2.16 ± 0.16 


0.129 



Table 2. Results of the fits to two resonating solutions of 5- wave intensity Is in 
7r~p — > TT~TT'^n at 17.2 GeV/c using the parametrization (4.1). 
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Is 






Ns 


xVd.o.f. 


Solution 


(MeV) 


(MeV) 






(1,1) 


732± 17 


244 ± 36 


0.60 ± 0.12 


0.913 


(1,2) 


700± 33 


300 ± 66 


1.21 ± 0.22 


0.126 


(2,1) 


740± 32 


296 ± 116 


1.02 ± 0.29 


0.204 


(2,2) 


711± 20 


300 ± 60 


1.70 ± 0.21 


0.238 



Table 3. Results of the fits to the four solutions of S-wave intensity Is in 7r+n 
TT+TT^p at 5.98 GeV/c using the parametrization (4.1). 
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Il 


rup 




r 


Nl 




Solution 


(MeV) 


(MeV) 


(GeV-i) 






(1,1) 


777± 1 


160 ± 2 


4.9 ± 0.5 


4.52 ± 0.03 


2.918 


(1,2) 


777± 1 


157 ± 3 


4.5 ± 0.7 


4.50 ± 0.05 


1.106 


(2,1) 


777± 1 


160 ± 3 


4.8 ± 0.7 


4.50 ± 0.04 


1.551 


(2,2) 


777± 1 


157 ± 3 


4.3 ± 0.6 


4.48 ± 0.04 


1.596 



Table 4. Results of the fits to the four solutions of P-wave intensity /l in tt p 
TT^TT+n at 17.2 GeV/c using the parametrization (4.1). 
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Il 


rup 




r 


Nl 




Solution 


(MeV) 


(MeV) 


(GeV-i) 






(1,1) 


780± 3 


195 ± 8 


5.9 ± 1.1 


1.64 ± 0.05 


0.670 


(1,2) 


779± 5 


189 ± 12 


4.3 ± 1.3 


1.59 ± 0.07 


0.214 


(2,1) 


778± 5 


192 ± 11 


4.7 ± 1.3 


1.59 ± 0.07 


0.206 


(2,2) 


779± 5 


182 ± 12 


3.6 ± 1.0 


1.55 ± 0.07 


0.212 



Table 5. Results of the fits to the four solutions of P-wave intensity II in 7r+n 
TT+TT^p at 5.98 GeV/c using the parametrization (4.1). 
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Figure Captions 

Fig. 1. Number of physical solutions out of 30 000 Monte Carlo selections of SDM 
elements as a function of dipion mass for iT~p — > 7r~n~^n at 17.2 GeV/c and 
-t =0.005-0.2 (GeV/c),2 and for ir+n Tr+ir-p at 5.98 GeV/c and -t =0.2- 
0.4 (GeV/c).2 

Fig. 2 Distributions of physical values of cos(7^^) for Tr+n — > tt+tt^p at 5.98 GeV/c 
and -t =0.2-0.4 (GeV/c)^ in the dipion mass bin 520 < m < 600 MeV. The 
distribution for Solution 1 yields a range of cos (7^^) from 0.52 to 1.00 with an 
average value of 0.94. The distribution for Solution 2 yields a range of cos(7g^) 
from 0.31 to 1.00 with an average of 0.64. 

Fig. 3. The mass dependence of physical solutions for moduli squared of iS-wave and P- 
wave nucleon transversity amplitudes and cosines of their relative phases in the 
reaction tt~p — > TT~n~^n at 17.29 GeV/c and momentum transfers —t =0.005- 
0.2 (GeV/c).^ The results are in the t-channel dipion helicity frame. 

Fig. 4. The mass dependence of physical solutions for moduli squared of 5'-wave and 
P-wave nucleon transversity amplitudes and cosines of their relative phases in 
the reaction Tr+n — > tt+tt^p at 5.98 GeV/c and momentum transfers —t =0.2- 
0.4 (GeV/c).^ The results are in the t-channel dipion helicity frame. 

Fig. 5. The mass dependence of physical solutions for moduli squared of S-wave 
transversity amplitudes S and S in n~p 7r~7r"'"n at 17.2 GeV/c and mo- 
mentum transfer —t =0.005-0.2 (GeV/c).^ Both solutions for the amplitude 
IS'pS resonate at 750-770 MeV while the amplitude l^pS is nonresonating in 
both solutions. 



Fig. 6. Four solutions for S'-wave partial wave intensity Is in the reaction tt p — > 
TT-TT+n at 17.2 GeV/c and -t =0.005-0.2 (GeV/c).^ The curves are fits to the 
Breit-Wigner form (4.1). The fitted parameters are given in Table 2. 

Fig. 7. Four solutions for ^-wave partial wave intensity Is in the reaction Tr+n 

TT+TT-p at 5.98 GeV/c and -t =0.2-0.4 (GeV/c).^ The curves are fits to the 
Breit-Wigner form (4.1). The fitted parameters are given in Table 3. 

Fig. 8. Four solutions for S'-wave partial wave intensity Is in the reaction 7r~^n — > 
TT+TT-p at 11.85 GeV/c and -t =0.2-0.4 (GeV/c).^ 

Fig. 9. Comparison of mass dependence of unnormalized moduli of P-wave production 
amplitudes and associated partial wave intensities in reaction n~^n^ — > n~^n~p 
at 5.98 GeV/c and -t =0.2-0.4 (GeV/c)? Shown are solutions 1. The other 
combinations of solutions are similar. The lines are to guide the eye. 
Fig. 10. Comparison of mass dependence of unnormalized moduli of P-wave production 
amplitudes and associated partial wave intensities in reaction 7r~p — > Tr^Tr+n 
at 17.2 GeV/c and -t =0.005-0.2 (GeV/c).^ Shown are solutions 1. The other 
combinations of solutions are similar. Also shown is the Briet-Wigner fit to 
spin-averaged partial wave intensity 1^(1,1). The Breit-Wigner fit is scaled 
and compared with the moduli of spin amplitudes \L\'^T, and |-ZvpS. 
Fig. 11. Breit-Wigner fit to spin-averaged partial wave intensity 1^(1,1) in Tr+n — > 
TT'TT+p at 5.98 GeV/c for -t = 0.2 - 0.4 (GeV/c).^ The Breit-Wigner fit is 
scaled and compared with the moduli of spin amplitudes and 
Fig. 12. Test of predictions p^g + Pqq -|- 2p^j^ = — 2(pqq — p\-^) = -\-2p\_^ due to vanishing 
^i-exchange in the reaction tt~p — > Ti~Ti^n at 17.2 GeV/c and —t =0.005-0.2 
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(GeV/c).2 

Fig. 13. Test of predictions Repfg = Rep\g = Rep^g — due to vanishing ^di-exchange 
in the reaction ir'p tt-tt+ti at 17.2 GeV/c and -t =0.005-0.2 (GeV/c).^ 

Fig. 14. ^-wave intensity normahzed to one at maximum value. The data correspond 
to solutions /s(l, 1) and /s(2, 2) at 17.2 GeV/c. The smooth curves are predic- 
tions based on tt+tt" — > tt+tt" analyses (A and D from Ref. 27, B from Ref. 32 
and C from Ref. 33) and on current algebra and PC AC (E from Ref. 34). The 
curves are taken from Ref. 15. Notice that the curves correspond to / = S'- 
wave intensity in tt+tt" tt+tt" reaction while the data correspond to ^-wave 
intensity in 7r~p — > 7r~7r+n in the physical region of t and include contribution 
from 1 = 2 S -waves. 

Fig. 15. The 1 = 2 contribution to S'-wave intensity in it^tt~ scattering. The solid 
curve is / = S-wave intensity in tt+tt" — > tt+tt" from Ref. 35. The dashed 
curve shows this ^-wave intensity with 1=2 contribution included. The 1 = 2 
phase shifts were taken from Ref. 37. The data correspond to solution /^(l, 1) 
and /s(2,2) at 17.2 GeV/c. 
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